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ABSTRACT 

We have extended our analytical chemical evolution modelling ideas for the Galaxy 
(Pagel & Tautvaisiene 1995, 1997) to the Magellanic Clouds. Unlike previous authors 
(Russell & Dopita 1992; Tsujimoto et al. 1995; Pilyugin 1996), we assume neither 
a steepened IMF nor selective galactic winds, since among the a-particle elements 
only oxygen shows a large deficit relative to iron and a similar deficit is also found 
in Galactic supergiants. Thus we assume yields and time delays identical to those 
that we previously assumed for the solar neighbourhood. We include inflow and non- 
selective galactic winds and consider both smooth and bursting star formation rates, 
the latter giving a better fit to the age-metallicity relations. We predict essentially 
solar abundance ratios for primary elements and these seem to fit most of the data 
within their substantial scatter. Our LMC model also gives a remarkably good fit to 
the anomalous Galactic halo stars discovered by Nissen & Schuster (1997). 

Our models predict current ratios of SNIa to core-collapse supernova rates en- 
hanced by 50 per cent and 25 per cent respectively relative to the solar neighbourhood, 
in fair agreement with ratios found by Cappellaro et al. (1993) for Sdm-Im relative to 
Sbc galaxies, but these ratios are sensitive to detailed assumptions about the bursts 
and a still higher enhancement in the LMC has been deduced from X-ray studies of 
remnants by Hughes et al. (1995). The corresponding ratios integrated over time up 
to the present are slightly below 1, but they exceed 1 if one compares the Clouds with 
the Galaxy at times when it had the corresponding metallicities. 

Key words: Magellanic Clouds; stars: abundances; stars: mass function; supernovae: 
general 



1 INTRODUCTION 

The histories of star formation and chemical evolution in 
the Magellanic Clouds exhibit distinct features in compari- 
son to those of the disk of the Milky Way. In both Clouds, as 
reviewed by Olszewski, Suntzeff & Mateo (1996) and West- 
erlund (1997), it seems that the majority of the stars are un- 
der 4 Gyr old, with a sprinkling of older stars. In the LMC 
there is a sudden rise in the star formation rate (SFR) 2 to 
4 Gyr ago (Elson, Gilmore & Santiago 1997; Geisler et al. 
1997) preceded by either a constant lower SFR (Geha et al. 
1997) or possibly a virtual gap as manifested by the cluster 
age distribution (Da Costa 1991; van den Bergh 1991), itself 
preceded by the formation of a minority older population 
resembling the Galactic globular clusters although having 
disk-like kinematics, while in the SMC the SFR appears to 
have been more uniform (van den Bergh 1991; Olszewski, 
Suntzeff & Mateo 1996). In the LMC, the sudden rise in the 
SFR a few Gyr ago is reflected in a corresponding sudden 



rise in the metallicity, both from [Fe/H] in the clusters (Ol- 
szewski et al.; Geisler et al.) and from a-particle elements 
in planetary nebulae (Dopita et al. 1997), while in the SMC 
there are fewer data, but some signs of a break in the age- 
metallicity relation (AMR) at an age between 2 and 4 Gyr 
(see below). From Olszewski et al., the relative numbers of 
known clusters older and younger than 3 Gyr are 15:100 in 
the LMC and 5:37 in the SMC, i.e. about equal proportions, 
but only one cluster is known in the LMC with age between 
3 and 10 Gyr (and that one could be an interloper), while 
in the SMC there are 3 clusters known in the same interval. 

Many models have been put forward for the chemical 
evolution of the Magellanic Clouds, paying attention to the 
distinct Fe/O and Fe/a ratios which are generally found 
to be higher than in Galactic stars with the same metallic- 
ity, i.e. Fe/H. Gilmore & Wyse (1991) pointed out that one 
way to get this effect is to assume separated star formation 
bursts, with SNIa contributing extra iron during quiescent 
intervals. Russell & Dopita (1992) put forward a very de- 
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tailed and comprehensive pair of models covering several 
element: element ratios previously determined by themselves 
and others and normalised to the solar neighbourhood (ISM 
and F supergiants) rather than the Sun. This model assumed 
a smooth evolution in which the Clouds are built up by in- 
flow of unprocessed gas according to certain time scales and 
a slightly steeper IMF than for the solar neighbourhood, 
which helps to bring down O/Fe and a/Fe. They noted that 
the observations seem to indicate [O/Fe] slightly lower than 
[a/Fe], which might be a sign of O/a stellar yield ratios 
increasing with mass. They found but minor variations in 
heavy-element:iron ratios relative to the solar neighbour- 
hood in the LMC, with the exception of an excess of Nd 
and Sm which have a significant contribution from the r- 
process. In the SMC they found a still larger excess of Nd 
and the essentially pure r-process elements Sm and Eu, av- 
eraging 0.6 dex, and concluded that the main s-process had 
been less effective, and the r-process more effective, than in 
the solar neighbourhood. Russell & Dopita's work referred to 
young objects (supergiants and the interstellar medium) and 
therefore did not include discussion of the age-metallicity re- 
lation. 

Another similarly comprehensive set of models was put 
forward by Tsujimoto et al. (1995), with the express purpose 
of explaining a relatively high SNIa/SNII ratio, deduced 
from the Fe/O ratio and supported by some direct counts of 
remnants, as had been proposed earlier by Barbuy, de Fre- 
itas Pacheco & Castro (1994). Tsujimoto et al. considered 
both smooth and bursting models. From an optimal fit to 14 
elemental abundances from O up to Ni and smooth (burst- 
ing) chemical evolution models they deduced values of 0.15, 
0.24 (0.21) and 0.3 (0.28) for the ratio, integrated over time, 
of SNIa/SN(II+Ib+Ic) in the solar neighbourhood, the LMC 
and the SMC respectively, and suggested that they origi- 
nate from Salpeter-like power-law exponents in the IMF of 
respectively 1.33, 1.71 (1.62) and 1.88 (1.84). These steeper 
IMFs (between fixed upper and lower mass limits) , also pos- 
tulated by Russell & Dopita, then also serve to produce the 
subsolar metallicities in the Magellanic Clouds at the cor- 
responding gas fractions and give reasonable age-metallicity 
relations, especially with the burst models. 

The hypothesis of a slightly steeper IMF compared to 
the Galaxy is neither supported nor ruled out by direct 
star counts (Hill, Madore & Freedman 1994 and references 
therein). An alternative way to explain low metallicities is 
to assume outflow (e.g. Carigi et al. 1995). Such outflow 
can be assumed to be either homogeneous or selective, the 
latter case being associated with starbursts and leading to 
enhancement of the Fe/O ratio, among others (e.g. Mar- 
coni, Matteucci &: Tosi 1994). Such a model has been put 
forward for the LMC by Pilyugin (1996), who points out 
that a closed bursting model, like those of Gilmore & Wyse, 
will produce the same final oxygen abundance at a given 
gas fraction as a smooth model, so that some other factor is 
needed to provide the reduction in metallicity compared to 
the solar neighbourhood with the same gas fraction. Given 
the absence of compelling evidence for a change in the IMF, 
Pilyugin postulates selective winds associated with star for- 
mation bursts; an initial square-wave burst in the first 0.2 
Gyr, accounting for 8 per cent of the stars, is followed by a 
hiatus until 9 Gyr, when a major burst or series of bursts 
begins, ending after 12.5 Gyr (0.5 Gyr before present) and 



accounting for the remainder, in accordance with the age 
distribution of the LMC clusters. A model with both ho- 
mogeneous and selective galactic winds (with a selective ex- 
pulsion of 75 per cent of ejecta from massive stars and a 
homogeneous expulsion of a mass equal to that being made 
into stars) is found to give good fits to the AMR and to 
[Fe/O] . 

Unfortunately, the assessment of [Fe/O] and, to a lesser 
extent, other abundance ratios is dependent on which Galac- 
tic standards are used (cf Pagel 1992). While it is certainly 
true that [Fe/O] is not negative in the Clouds, it is not 
strongly positive either when compared to Galactic super- 
giants instead of the Sun, and because of various system- 
atic effects the supergiants may provide a better standard. 
Among a-particle elements, notably Mg and Ca, and for Ti 
which seems to behave like the a-elements among Galac- 
tic stars, the offset from either solar or local supergiant ra- 
tios is insignificantly different from zero, compared to small 
positive values in the Milky Way for the same [Fe/H] (Ed- 
vardsson et al. 1993). There is thus some difference, but not 
necessarily large enough to demand either a changed IMF 
or selective galactic winds. Some sort of winds are prob- 
ably present, as witness the Magellanic stream, and these 
may well account for the relatively low metallicities, and we 
agree with Russell & Dopita and with Tsujimoto et al. in 
assuming a significant effect of inflow, since this is needed, at 
least in smooth models with reasonable star formation laws, 
to account for the 'G-dwarf problem posed by the relatively 
small number of old objects. The existence of separated star 
formation bursts is well attested in the LMC by population 
studies of both clusters and field stars, and is supported in 
both Clouds by the AMRs. 

With this background, it seems worth while to inves- 
tigate how well the chemical evolution of the Magellanic 
Clouds can be modelled with yields and time-delays equal 
to those applying to the solar neighbourhood, but with in- 
flow and a homogeneous wind. In two previous papers (Pagel 
& Tautvaisiene 1995, 1997) we developed analytical models 
giving a fair fit to stellar abundances of primary elements in 
the solar neighbourhood, and in this paper we shall attempt 
to apply the same ideas to the Clouds. 

2 THE MODEL 

2.1 General description 

In agreement with Russell & Dopita (1992) and Tsujimoto 
et al. (1995), we assume the Clouds to have been built up by 
gradual infall of unprocessed material. This helps to alleviate 
the 'G-dwarf problem that would otherwise arise from the 
small relative proportion of old, metal-poor stars. Our for- 
malism is an adaptation of that of Pagel (1997) as detailed 
below; like Tsujimoto et al, we assume linear laws of star 
formation and investigate both smooth and bursting mod- 
els. Like Pilyugin (1996), we assume yields and time delays 
identical to those which apply to the solar neighbourhood 
and appeal to galactic winds to explain the low metallicities 
of the Clouds in relation to their current gas fractions; we 
take the specific numbers from our two previous papers. Un- 
like Pilyugin, however, we ignore selective winds and assume 
just a non-selective wind proportional to the SFR, similar 
to the model of Hartwick (1976) for the Galactic halo. 
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Figure 1. SFR history for the LMC according to our model. The 
full curve shows the bursting model according to Table 1, while 
the broken-line curve shows a smooth model with u = 0.18t. 
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Figure 2. SFR history for the SMC according to our model. The 
full curve shows the bursting model according to Table 2, while 
the broken-line curve shows a smooth model with u = 0.115t. 



2.2 The formalism 

We assume a linear or quasi-linear star formation law 

ds _ 
du ^' 

where s is the mass in stars (+ remnants), g is the mass 
of gas and u = J* uj(t')dt' where t is time and lo is the in 



(1) 



verse time scale for star formation, assumed constant in the 
'smooth' models. In 'bursting' models, u) is assumed to be 
constant over certain time periods, between which it changes 
discontinuously. Inflow is assumed to occur at a rate 



/(t)= W (t)e-», 



(2) 
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Figure 3. Age-abundance relation for a-elements in the LMC. 
Curves are as in Fig. 1. Data points are from Dopita et al. (1997): 
open circles represent oxygen while filled circles represent an av- 
erage of Ne, S and Ar. 



and outflow at a rate 
ds 

6 ^ = 11 dt = VUJ9, ^ 
with T) = const., so that the gas mass and total mass satisfy 
the differential equations 

dg 
du 
and 
dm 



(4) 



du 



V9 



(5) 



respectively. From Eqs. (4), (5), the gas and total masses 
evolve according to 

g(u)= [e- u ~e-^]/ V ; (6) 
m{u)= [l-e-d+^j /(! + „); (7) 
the star mass is 



s(u) = m(u) — g(u) 
and the gas fraction is 
jx(u) = g(u)/m(u). 



(8) 



(9) 



Using the instantaneous recycling approximation, the abun- 
dance of a promptly produced primary element such as O 
or Mg, in units of its yield (assumed constant) , satisfies 



^(321) = - (1+77)21], 
which has the solution 



(10) 



(11) 



zi(u) = - 
V 

Results of Eqs. (1) to (11) for our assumed model para- 
meters are given in the first seven columns of Tables 1 and 
2. Figs 1 and 2 show the resulting SFR histories for both 
smooth and bursting models; the rough coincidence in time 
between our assumed bursts in the two Clouds suggests that 
they may have resulted from some mutual interaction. For 
the LMC, our assumed value of 1 for rj agrees with that 
adopted by Pilyugin and the current star formation rate is 
about average; this may be compared with the remark by 
Westerlund (1997) that the current SFRs in both Clouds 
deduced from Ha emission may indicate that their current 
SFRs are below their averages. Thus the decline that we 
assume in the last 3 Gyrs from the peak of the recent burst, 
resulting from our assumption of Eq. (1) with a constant u 
in the relevant interval, seems to be fairly realistic. The final 
gas fraction adopted for the SMC is in accordance with the 
number given by Westerlund; for the LMC it is somewhat 
larger than his figure of < 8%, but in accordance with that 
adopted in other models in the literature and similar to some 
estimates for the solar neighbourhood. Fig 3 shows the age- 
abundance relation in the LMC for O, Ne, S and Ar, assumed 
to be instantaneously produced with a yield of 0.7 times 
their solar abundance, and confirms that a bursting model 
gives a better fit than a smooth one. 

For elements (or components thereof) that are ejected 
after a significant time delay, we use the 'delayed production' 
approximation (Pagel 1989) which assumes the element to 
be ejected at a fixed time delay A after the time of star 
formation, whereas the total mass ejection from a generation 
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Table 1. Bursting model for the LMC; (u>) = 0.18 Gyr ~ 1 ;n = l. 

u> t u g m fi zi 22 22 22 22 [O/H] [Fc/H] 

Gyr" 1 Gyr 0.023 0.037 1.33 2.67 
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0.003 
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00 
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0.000 
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-2 


.71 


-3.11 


0.500 


0.075 


0.067 





.070 


0. 


96 





.037 


0. 


.034 


0.032 


0.000 


0.000 


-1 


.59 


-1.98 


0.900 


0.135 


0.110 





.118 


0. 


93 


0. 


.066 


0. 


063 


0.061 


0.000 


0.000 


-1 


.34 


-1.73 


1.330 


0.200 


0.148 





165 





90 


0. 


096 


0. 


.094 


0.092 


0.000 


0.000 


-1 


.17 


-1.57 


2.000 


0.300 


0.192 





.226 


0. 


85 





.143 


0. 


140 


0.139 


0.022 


0.000 


-1. 


.00 


-1.31 


2.023 


0.302 


0.193 





.227 





85 


0. 


.143 





.143 


0.141 


0.024 


0.000 


-1 


.00 
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2.037 


0.303 
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85 





.144 





.143 


0.143 


0.025 


0.000 


-1 


.00 


-1.30 


2.670 


0.354 


0.209 





.254 


0. 


83 





.166 





.166 


0.165 


0.073 


0.000 


-0 


.93 
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3.330 


0.406 


0.222 





.278 


0. 


80 





.190 


0. 


.189 


0.188 


0.135 


0.007 


-0 
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4.670 


0.514 
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.321 





75 


0. 


.235 


0. 
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0.046 
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-0.84 


6.000 
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70 
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.377 
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.395 
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63 





.340 





.339 
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0.306 


0.215 


-0 


.62 


-0.65 


9.000 


0.860 


0.244 





411 
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.369 





.369 
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-0.61 


10.000 
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-0 


56 
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.435 
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0.349 
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1.032 
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11.037 
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-0 
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1.685 
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.27 


0. 


007 


0. 


657 


0.657 


0.455 


0.234 


-0 


.35 


-0.43 


13.000 


2.020 


0.115 





491 


0. 


23 


0. 


691 


0. 


693 


0.695 


0.626 


0.239 


-0 


.32 


-0.34 


13.670 


2.355 


0.086 





.496 


0. 


.17 





.753 


0. 


.757 


0.760 


0.915 


0.270 


-0 


.28 


-0.23 


13.800 


2.420 


0.081 





.496 


0. 


16 





.764 





.769 


0.771 


0.963 


0.422 


-0 


.27 


-0.21 


14.000 


2.520 


0.074 





.497 





15 


0. 


780 





.785 


0.787 


1.033 


0.645 


-0 


.26 


-0.19 



0.50 

Table 2. Bursting model for the SMC; (w) = 0.115 Gyr" 1 ; r) = 2. 

u> t u g m fi 2i 22 22 22 22 [O/H] [Fc/H] 

Gyr" 1 Gyr 0.023 0.037 1.33 2.67 
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Figure 4. Age-mctallicity relation for the LMC. Curves represent 
our models as in previous figures. Data sources are as follows: open 
circles, Olszewski et al. (1991); crosses, Girardi et al. (1995); open 
squares, Geisler et al. (1997); filled circles, Idiart & de Frcitas 
Pacheco (1997). 




Figure 5. Age-metallicity relation for the SMC. Curves represent 
our models as in previous figures. Data represented by open circles 
are from Da Costa (1991). 
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Thus in separate time segments A to ti (or t\ to A if 
A > ti), ti (or A) to h+A, ti+A to ti, t 2 to t 2 + A, t 2 + A 
to £3 —T,u)' — const, and the exponents in Eq. (13) are all 
linear functions of time: 



(1 + 77)11 
u 



p + qt; 
r + st, 



(20) 



say, and so Eq. (13) can be solved analytically for each in- 
terval t a to £(,: 
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P-{1 + V)r 

- (l + »7> 
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of stars is still assumed to take place instantaneously. In this 
case the abundance 22 in units of the yield varies according 
to 

d 



{922 



1 f 

■■ w g 



(1 + rj)z 2 ujg; t > A, 



(12) 



dt 

where the dashes indicate that the variable has to be taken 
at the time t — A rather than t. Eq. (12) leads after some 
reduction to 



|[»(e--l)] = 

tion 

zi{u) 



-(1 + Tj)u'l 



(13) 



For smooth models (u> — const.), Eq. (13) has the solu- 

e i(«-» 4 ) _ 1 



e »7« _ 1 

U > UjA. 



V 



- (u - wA) 



(14) 



In bursting models, we assume ui = uii , say, up to time 
ti, UJ = 0J1 between t\ and and lo — lo^ between £2 and 
£3 = T, the age of the system, assumed to be 14 Gyr. L02 
is small, representing a quiescent phase between an initial 
starburst represented by ui and a stronger recent starburst 
represented by 0*3. Then, in Eq. (13), 



J = u = 0; for t < A; 

J = wi; 

u — uj\(t — A); for 



(15) 



2.3 Yields and time delays 

Table 3 gives the yields and time delays that we assume, 
after Pagel & Tautvaisiene (1995, 1997). Yields are in units 
of solar abundance of the corresponding element and time 
delays in Gyr are given in the top row. Total abundances 
are obtained for each time step by summing the products of 
these yields with the appropriate values of z for the relevant 
time delay, which are given in columns 7 to 11 of Tables 1 
and 2. Figs 4 and 5 show the resulting age-metallicity rela- 
tions for the two Clouds, which are similar to those derived 
by Tsujimoto et al. (1995), but differ somewhat from those 
for the LMC by Pilyugin (1996). 



3 ELEMENTiELEMENT RATIOS 

The last seven columns of Tables 1 and 2 give abundances 
or constituent z- values for individual elements in the Clouds 
calculated from our bursting models, and elementdron ra- 
tios are plotted in Figs 6 to 11. The observational data, 
which come mainly from MC supergiants, have in some cases 
been adjusted to be relative to Galactic supergiants rather 
than the Sun, since (as noted by Russell & Dopita 1992, 
Pagel 1992 and Hill, Andrievsky & Spite 1995) there are 
some discrepancies, particularly in the case of oxygen, and 
it seems better to compare like with like. Russell & Dopita 
took [O/Fe] for the solar neighbourhood to be —0.03 only, 
but various studies of Galactic supergiants (Luck & Lambert 
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Figure 6. Elcment:iron ratios for oxygen, a-particle elements and 
Ti in the LMC and in the anomalous Galactic halo stars. The full- 
drawn curves show predictions from the bursting model and the 
broken-line curves those from the smooth model. Data points for 
the LMC stars have been shifted upwards by 0.2 dcx for oxygen 
and downwards by 0.1 dcx for silicon. Data sources: + signs, Luck 
& Lambert (1992); open circles, Thevcnin (1997); five-cornered 
stars, Hill, Andriewski & Spite (1995); open squares, Spite, Bar- 
buy & Spite (1993); open diamonds, Jiittncr et al. (1992); open 
triangles, McWilliam & Williams (1991); small six-pointed stars, 
Barbuy, de Freitas Pacheco & Castro (1994); asterisks, Richtler, 
Spite & Spite (1989); crosses, Russell & Bessell (1989); large six- 
pointed stars, anomalous Galactic halo stars after Nissen & Schus- 
ter (1997). 
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Figure 7. Element:iron ratios for s-process elements in the LMC 
and in the anomalous Galactic halo stars. Curves and data sources 
as in Fig 6. Data points for Ba and La in the LMC stars have 
been shifted downwards by 0.1 dex and 0.16 dex respectively. 

1985; Spite, Barbuy & Spite 1989; Hill, Andrievsky & Spite 
1995; Venn 1995) indicate a more negative value, [O/Fe] 
~ —0.2 relative to the Sun. The observational [O/Fe] ra- 
tios plotted in Figs 6 and 9 have accordingly been shifted 
upwards by 0.2 dex, which leads to a somewhat different pic- 
ture from that assumed in previous chemical evolution mod- 
els for the Clouds; in particular, there is now (apparently) 
good agreement between oxygen and the other Q-elements. 
The oxygen abundances from planetary nebulae shown in 
Fig 2 also support this upward adjustment. On the basis 
of studies of Canopus and other Galactic supergiants by 
Luck (1982), Reynolds, Hearnshaw & Cottrell (1988), Rus- 
sell & Bessell (1989) and Spite, Spite & Frangois (1989), wc 
have shifted the data points for [Si/Fe], [Ba/Fe] and [Nd/Fe] 
downwards by 0.1 dex and [La/Fe] downwards by 0.16 dex, 
taking all the remaining elements at face value. 

Although the model curves in Figs 6 to 11 are meant 
to represent an evolution in time, it has to be borne in 
mind that the data points come from young supergiants, so 
that their horizontal spread is due to scatter (part of which 
is probably observational) and not from evolution in time. 
Consequently it is only the centroid of these points that is 
significant and no evolutionary trend can be deduced from 
the data. However, Nissen & Schuster (1997) have found a 
sample of stars from the outer halo of our Galaxy in which, 
among other peculiarities, the a-element to iron ratios are 
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Figure 8. Elemcntiiron ratios for heavy metals in the LMC. 
Curves and data sources as in Fig 6. Data points for Nd have 
been shifted downwards by 0.1 dex. 

more or less solar rather than enhanced as in other stars with 
the same metal deficiency. They suggest that these stars 
could have been captured from dwarf galaxies such as the 
Magellanic Clouds and we accordingly include them (with- 
out any adjustments) in the relevant plots. Our LMC model 
fits these ratios very well, including a slight rise from solar- 
like [q/Fc] around [Fe/H] = -0.8 to [a/Fe] ~ 0.2 at [Fe/H] 
= —1.2 (Fig 6) and a more uniform trend for s-process el- 
ements (Fig 7), supporting an origin in the LMC (or some 
similar system). The fit of our SMC model (Figs 9, 10) is 
not quite as good. As far as the actual Magellanic stars are 
concerned, there is considerable scatter in the data, much of 
which we suspect to be unreal, although we have chosen to 
plot individual determinations rather than to follow Russell 
& Dopita in just plotting an average. Within the scatter, 
the agreement of our model with abundances deduced from 
observation seems to be quite satisfactory, with the pos- 
sible exception of a few of the heaviest elements. Broadly 
speaking, there is hardly any evidence for significant depar- 
tures from solar, or solar-neighbourhood, abundance ratios 
among the elements considered (we do not discuss carbon, 
nitrogen or sodium which all exhibit 'secondary' behaviour), 
and that result is just what is predicted; whereas the steep- 
ened IMF models of Tsujimoto et al, while successful for 
oxygen (if unadjusted to Galactic supergiants) , predict sub- 
solar abundances (relative to iron) of several a-particle el- 
ements, notably magnesium, which are not supported by 




-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

[Fe/H] 



Figure 9. Elcmcntiiron ratios for oxygen, ct-particlc elements 
and titanium for the SMC and anomalous Galactic halo stars. 
Data points for oxygen and silicon in the SMC have been shifted 
upwards by respectively 0.2 dex and 0.1 dex . Solid and broken 
curves represent our SMC bursting and smooth models respec- 
tively. Data sources: open diamonds, Jiittner et al. (1992); filled 
diamonds, Hill, Barbuy & Spite (1997); squares with + signs, 
Spite, Barbuy & Spite (1989); open squares, Hill & Spite (1997); 
open triangles, Spite et al. (1986); open circles, Spite, Richtlcr 
& Spite (1991); filled circles, Thevenin (1997); + signs, Luck 
& Lambert (1992); crosses, Russell & Bessell (1989); asterisks, 
Spite, Spite & Francois (1989); filled triangles, Hill (1997); filled 
five-pointed stars, Thevenin & Foy (1986); open five-pointed stars, 
Foy (1981); open six-pointed stars, anomalous Galactic halo stars 
from Nisscn & Schuster (1997). 

the data. Pending improved data, especially for the heaviest 
elements, we conclude that neither a steepened IMF nor se- 
lective galactic winds are required to explain the abundances 
in the Magellanic Clouds. 



4 RELATIVE SUPERNOVA RATES 

Barbuy, de Freitas Pacheco & Castro (1994) and Tsujimoto 
et al. (1995) have claimed that the overall ratios of Type 
la supernovae to core-collapse supernovae ever formed in 
the Clouds are higher than in the solar neighbourhood. Our 
models do not support this and they fall down if it is true. 
With our approximations, the ratio is proportional to 

s(T-A) m(T - A) - g(T - A) 

s(T) m(T)-g(T) ' 1 > 
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Figure 10. Elcmcnt:iron ratios for s-proccss elements in the SMC 
and anomalous Galactic halo stars. Curves and symbols as in Fig 
9. Data points for Ba and La in the SMC stars have been shifted 
downwards by 0.1 dex and 0.16 dex respectively. 
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Figure 11. Element:iron ratios for heavy metals in the SMC. 
Curves and symbols as in Fig 9. Data points for Nd have been 
shifted downwards by 0.1 dex. 



Taking u = 4.5, ujA — 0.4 for the solar neighbourhood 
(Pagel & Tautvaisiene 1995), this ratio comes out to be 0.95, 
whereas in the LMC and SMC according to Tables 1 and 2 
it is 0.86 and 0.77 respectively, i.e. marginally lower. The 
situation changes if we compare the above Magellanic ra- 
tios with those that prevailed in the solar neighbourhood at 
those times in the past when it had the corresponding metal- 
licities; these turn out to be 0.81 for [Fe/H] = -0.2 and 0.62 
for [Fe/H] = —0.45. Thus we attribute the lower a/Fe ratios 
chiefly to the effect of slower evolution up to a metallicity 
which prevailed in our Galaxy at an earlier time when there 
had indeed been relatively fewer SNIa in the latter. It can 
be seen from Figs 6 to 11 that the influence of starbursts on 
the element:element ratios is comparatively insignificant in 
our models. 

When we consider current supernova rates, the situation 
is quite different again. In this case the ratio is proportional 
to 

(ds/dt) T -A g(T - A) 
(ds/dt) T g(T) ■ ( > 

The corresponding numbers are 1.2 for the solar neighbour- 
hood, 1.8 for the LMC and 1.5 for the SMC in our bursting 
models. Thus we expect a 50 per cent higher ratio of SNIa to 
core-collapse supernovae in the LMC at present. Assuming 
a ratio in the solar neighbourhood between 0.11 and 0.25 
(van den Bergh & McClure 1994), the ratio in the LMC is 
then predicted to be between 0.16 and 0.38, which is in good 
agreement with the average figure of 0.29 given for Sdm-Im 
galaxies by Cappellaro et al. (1993). A still higher ratio, of 
order 1, has been estimated for the LMC on the basis of X- 
ray observations of remnants by Hughes et al. (1995), who 
suggest a likely lower limit of 0.25, so that our model agrees 
qualitatively, if not quantitatively, with their findings. Cur- 
rent ratios are sensitive to the precise assumptions made 
about the bursts, and do not necessarily represent historic 
averages. We suggest that the relatively high ratio of SNIa 
to core collapse SN observed in the LMC and other Sd-Im 
galaxies is related to their star formation history rather than 
the IMF. 
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